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ABSTRACT 

We show that gaseous Hi disks of primordial composition irradiated by an 
external radiation field can develop a multiphase medium with temperatures 
between 10 2 and 10 4 K due to the formation of molecular hydrogen. For a given 
Hi column density there is a critical value of the radiation field below which only 
the cold Hi phase can exist. Due to a time decreasing quasar background, the 
gas starts cooling slowly after recombination until the lowest stable temperature 
in the warm phase is reached at a critical redshift z = z CT . Below this redshift 
the formation of molecular hydrogen promotes a rapid transition towards the 
cold Hi phase. We find that disks of protogalaxies with 10 20 £ N HI < 10 21 
cm -2 are gravitationally stable at T ~ 10 4 K and can start their star formation 
history only at 2, after the gas in the central portion of the disk has 

cooled to temperatures T & 300 K. Such a delayed starbust phase in galaxies of 
low gas surface density and low dynamical mass can disrupt the disks and cause 
them to fade away. These objects could contribute significantly to the faint blue 
galaxy population. 

Subject headings: Galaxies: evolution, starburst - Molecular processes - 
Instabilities 
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1. Introduction 

One of the most interesting and controversial discoveries in the field of galaxy formation 
and evolution is related to the excess number count of Faint Blue Objects (see for example 
Koo 1996 or Ellis 1997 for a review). Various evolutionary models favour a scenario in which 
massive galaxies form stars at high redshifts (z ~ 2-3), while low-mass systems (dwarfs) 
experience an initial and disruptive starburst much later and contribute significantly to the 
Faint Blue Objects counts (Corbelli & Salpeter 1995, hereafter CS; Ferguson & McGaug 
1995; Babul & Ferguson 1996; Gwyn & Hartwick 1996; Cowie et al. 1997; Guzman et al. 
1997). The delayed star formation in dwarf galaxies has been explained as a consequence 
of the delayed recombination of ionized hydrogen (Babul & Rees 1992; Efstathiou 1992) 
and of the late onset of the gravitational instability which, in low mass systems, takes place 
only when the temperature drops well below 10 4 K (Babul & Rees 1993, CS). 

Cooling below 10 4 K can be caused by a slow build-up of metals, due to some sporadic 
star formation events before the onset of the large scale starburst (CS). However, since 
the metal abundance required for efficient cooling increases as the Hi column density of 
protodisks decreases, it is likely that additional processes must regulate the cooling in 
systems which are at lower end of the gas surface density distribution. In this Letter we 
examine if the formation of molecular hydrogen can provide such cooling in gaseous disks. 
The role of molecular coolants has been investigated in various cosmological contexts, such 
as the collapse and fragmentation of Jeans unstable clouds, pregalactic shocks and galaxy 
formation models (Shapiro & Kang 1987; Vietri & Pesce 1995; Ostriker & Gnedin 1996; 
Haiman et al. 1995, 1996; Anninos & Norman 1996; Tegmark et al. 1997; Kepner et 
al. 1997). All these studies have emphasized the importance of molecular hydrogen as the 
major agent to drive the temperature of the gas well below 10 4 K. The ability of H 2 to 
form in sufficient amount depends critically on the shielding of the gas to the background 
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ionizing and dissociating radiation. Here, we study the evolution of the central regions 
of a protogalactic disk which at some redshift are warm and in thermal equilibrium with 
an external radiation field. We show that, as the radiation field decreases with time, 
the formation of a small concentration of H 2 drives the gas out of equilibrium causing a 
rapid transition of the gas from the warm to the cold phase at a critical redshift z CT . The 
consequent reduction of thermal support can cause a gravitational instability of the disk 
and we examine which are the most favorable conditions for this to happen. 

2. Physical Background 

We consider low gas surface density proto-disks (LSPs) embedded in a spherical dark 
matter halo and rotationally supported against gravitational collapse. Similarly to what is 
today the disk of a low-surface brightness galaxy (de Blok et al. 1996), LSPs are expected 
to have slowly rising rotation curves, with lower asymptotic values, V^,, and shallower gas 
radial distributions compared to disks with higher gas surface density. We concentrate 
on the central regions of LSPs at it! £ R , with R in the range 3-8 kpc, where we can 
consider the dark matter density and the gas column density as radially constant. Here the 
column density of neutral hydrogen (which is by far the dominant species) is typically in 
the range 10 20 £ N HI <, 10 21 cm -2 , and the corresponding gas masses inside R are between 
10 7 -10 9 M Q . The dark matter determines the value of V^,; for = 30-150 km s _1 the 
total mass inside Ro is 10 8 -10 10 M . Disks with such characteristics are expected to exist 
already at z ~ 2, as shown by the theoretical calculations by Kauffmann (1996) and by 
observations of damped Lya systems (Wolfe et al. 1994). 

After recombination the vertical gas distribution is nearly isothermal, in hydrostatic 
equilibrium at each radius. It will be shown in the next Section that the gas stays quite 
uniform in temperature for most of the time and therefore in this Letter we do not consider 
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explicitly its vertical stratification. We approximate the total gas pressure with its value at 
midplane as 

Z ~ 10 f Nm ) 2 + 2 ( Nm \ ( Cs ) ( V °° ) (^Y 1 cm- 3 K (2 1) 
k ~ 40 20 cm- 2 ^ + Wo^cm-^Vkms-^kms^kpcJ 1 ' 

where k is the Boltzmann constant, and c s is the isothermal sound speed. The two 

terms in eq. (2.1) represent the contributions of the gas self-gravity and the dark matter, 

respectively. The pressure is radially constant in the central region and changes with time 

as the slab cools due to the temperature dependence of c s , which we follow explicitly. 

The chemistry of H 2 depends upon the density of free electrons and therefore on the 
ionizing flux at energies E > 100 eV, able to penetrate Hi column densities £ 10 20 cm -2 , as 
well as on the UV flux below the Lyman limit. At redshifts z <2, the quasar background 
decreases with time and accounts for most of the ionizing flux (Irwin, McMahon, & Hazard 
1991). We shall use the following fit to the results of Haardt & Madau (1996) at z — 1 

ixl.7xl0- 21 (j^)~°' 3 for E< 4.96 
A x 1.5 x 10- 22 (y^)~ 2 ' 4 for 4.96< E <13.6 



2 



13.6 



B x 5.0 x lQ-™(-£r) °' 3 for 13.6< E < 413 



S x 1.8 x 10- 24 (J-) L5 for£>413 
where E is in eV, and J(E,z) in erg s" 1 cm' 2 sr -1 Hz -1 . For the Haardt & Madau (1996) 
spectrum, A = B = 1 which defines our standard flux spectrum. We consider variations 
of A and B to account for uncertainties in the intrinsic emission spectrum of quasars, in 
the attenuation by intervening clouds, and in additional UV radiation from nearby bright 
star-forming galaxies or from sporadic internal star formation. 

We have followed the evolution of e, H, H + , H~ and H 2 by means of a reduced 
chemical network which includes the most relevant reactions listed in Table 1. We have 
also considered the H^ channel for H 2 formation, but its contribution was found to be 
negligible. Since temperatures of interest in this paper are below 10 4 K, secondary electrons 
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from H and He have been included following the prescription of Shull & van Steenberg 
(1985). These electrons tend to level the difference between He and H fractional ionization 
deriving from direct photoionization alone. Thus, the residual abundance of He + is set 
equal to that of H + . For the photoionization rates we have integrated the cross sections 
given in the references listed in Table 1 over the background spectrum. We account for the 
attenuation of the flux as it penetrates the slab to the midplane assuming that the density 
of the chemical species is uniform. The recombination coefficient for hydrogen includes 
all the recombinations to n > 2 (on the spot approximation) and we assume that all the 
recombinations of He produce one ionizing photon for H. The line-shielding of H 2 molecules 
in the ground level has been computed using the prescriptions of Federman et al. (1979) 
in the absence of dust grains. The photodissociation probabilities and oscillator strengths 
for Lyman and Werner bands are from Dalgarno & Stephens (1970) and Abgrall & Roueff 
(1989). 

The thermal evolution of the gas is computed by solving the energy equation together 
with the chemical network (we have used a Hubble constant of 75 km s" 1 Mpc" 1 ). Heating 
of the gas is provided by photoionization of atomic H and He (Shull & van Steenberg 1985). 
Cooling is due to electron-impact excitation of atomic H (Spitzer 1978) and to H-impact 
excitation of roto-vibrational levels of H 2 (Hollenbach & McKee 1989, but with the correct 
sign as in Flower et al. 1986). We have neglected the contribution of H 2 -H 2 collisions to the 
cooling because of the small molecular abundance. 

3. The two phase equilibrium and the fast gas cooling 

For a given value of V^/Rq the gas temperature of a slab in thermal equilibrium 
depends on the intensity of the background field and on the gas column density of the 
slab. For the standard background spectrum Fig. la (curves 1 and 2) shows the possible 
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equilibrium temperatures, after the gas has recombined, as function of redshift for two 
representative values of N HI . 

We see from Fig. la that conditions for a multiphase medium, in which gas at different 
temperatures coexists in pressure balance, can be achieved over a significant range of 
redshift. This is analogous to what happens in the present-day interstellar medium (Field 
et. al 1969), where metal lines radiation rather than H 2 dominates the cooling. In general 
a multiphase medium is possible if for a given pressure the cooling per atom is a non 
monotonic function of T. In the present context at T >, 7000 K the cooling rate is an 
increasing function of T since the dominant process is the e-H impact excitation. Below 
7000 K, H — H 2 impact excitation becomes the relevant cooling process: H 2 is formed 
efficiently by reaction (3) and is destroyed by reaction (5). The strong temperature 
dependence of reaction (5), proportional to e _2130 °/ T , allows a rapid rise of the H 2 fractional 
abundance as the temperature diminishes. Due to the weak temperature dependence of the 
H 2 cooling function in this region, cooling per particle becomes a decreasing function of 
T and the equilibrium curve has an inflection. This is a thermally unstable regime which 
ends at T £ 2000 K, when the 2-step photodissociation takes over reaction (5) (direct 
dissociation by hard photons, reaction 8, never contributes significantly to the H 2 balance). 
Under some circumstances it may not be obvious to determine the dominat thermal phase 
if, for a given ionizing flux, three equilibrium temperatures are possible. For the warm 
slab considered here the temperature evolution is much simpler: as the background flux 
declines with time, the temperature decreases monotonically from ~ 10 4 K moving along 
the equilibrium curve towards the right hand side of Fig., la. The non monotonic relation 
T(z) implies that at a certain redshift z cr the gas reaches the lowest stable equilibrium 
temperature in the warm phase: we shall call this transition point (marked by an asterisk 
in Fig. 1). For a slab of typical column density N HI = 10 20 4 cm" 2 z CT is ~ 0.5. The gas 
is thereafter forced to get out of equilibrium. The dashed tracks of Figure la show the 
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subsequent time evolution of T computed for a parcel of gas close to the midplane. We can 
see that the gas cools quickly toward the cold stable phase spending a very short interval 
of time at temperatures between 7000 K and 200 K. The transition towards the cold 
phase propagates rapidly vertically since the upper layers loose thermal support and settle 
towards the center incresing their density and cooling further. Due to the rapid cooling, 
shortly after the cold core forms most of the gas will be at close temperatures and the warm 
atmosphere left above contains little mass. During the transition the fractional abundance 
of H 2 , /(H 2 ), increases from ~3x 10~ 5 to ~ 10~ 3 and then levels off (Fig. lb). The final 
value of /(H 2 ) is not very sensitive to either the background flux or the gas column density. 

LSPs with N H i column density between 10 2(U and 10 208 cm -2 make a transition to 
the cold phase in the redshift interval z = 2 to z = 0. The value of z cr can vary appreciably 
with the intensity of the background flux. From curve 3 of Fig. 1(a) we can see in fact that 
as A increases, cooling of the slab is delayed until a much lower z. We have examined cases 
with 0.1 < A < 10 and 0.3 < B < 3 and the important point is that there is always a range 
of column densities between ~ 10 20 and ~ 10 21 cm -2 which make the transition to the 
cold phase in the redshift interval z = 0-2. This conclusion holds even if a small amount 
of metals (up to 0.01 solar) is present in the gas mixture as a contamination from an early 
generation of stars. Therefore, for the central regions of LSPs, molecular hydrogen plays a 
fundamental role in promoting and driving the thermal evolution of the gas. 



4. Gravitational instability and the star formation phase 

Protodisks exposed to a time-varying external radiation field follow the evolution 
described in the previous section if they are still gaseous and gravitationally stable in the 
warm phase. We examine the global gravitational stability of the simple hydrostatic model 
described in Sect. 2 using the Toomre criterion (Toomre 1964; see also CS). For R < Rq the 
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parameter Q can be written as 
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N H i 
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10 20 cm 



10 20 cm -2 
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(4.1) 



To evaluate the contribution of the self-gravity term we have used eq. (25) of Toomre (1963) 
assuming that the protogalaxy has a total gas mass M g ~ 4 x 10 8 (iViy//10 20 ) M Q . This 
guarantees a radially constant gas surface density for R < Rq and a slowly rising rotation 
curve with V^/Rq determined by the amount of dark matter. Therefore, for any N HI , we 
can derive Q as function of redshift. We are interested in all the disk models which are 
stable in the warm phase but become unstable (Q < 1) only at z < z CT due to the decrease 
of the isothermal sound speed driven by H 2 cooling. These galaxies will experience a delay 
in the onset of large scale star formation. 

A comprehensive view of the stability and dominant thermal phases of the gas in the 
plane N HI — V^/Rq for the standard background spectrum is given in Figure 2. In the 
region labelled cold & unstable slabs which are stable for z < z cr become gravitationally 
unstable before z — 0. LSPs in this region span a considerable range of Hi column densities 
below 10 21 cm' 2 . Models in the warm & unstable region have Q < 1 for T = 7000 K and 
the collapse starts as soon as the gas settles into the disk. The shaded regions in the left 
part of Fig. 2 show stable models which can be either cold or warm but cannot form stars. 
For some of these slabs cooling of the gas did actually start, but the temperature was not 
low enough to trigger the instability before z = (cold k, stable). Since cold Hi disks have 
not been observed, it is likely that LSPs were never formed with these initial conditions 
(very low N HI and high V^/Rq values). Finally, in the warm & stable region the gas 
remains warm up to z = 0; here the parameters of the slabs are more appropriate for the 
outer regions of disk galaxies (Corbelli & Salpeter 1993). Outer regions of LSPs and their 
central regions prior to z CT might also contribute significantly to QSO absorption lines (see 
also Linder 1997) because the time the gas spends in the warm phase at T £ 7000 K is long 
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compared to the transition time toward the cold and disruptive phase. 

The dotted lines in the cold & unstable region of Fig. 2 represent the loci of Q = 1 and 
are labelled with the corresponding redshift (the z = line coincides with the heavy solid 
line at the left hand side of this region). Systems of increasingly lower Hi column density 
start to cool and become unstable at lower and lower redshifts. For values of A and B 
higher than the standard case, the left border of the cold & unstable region and the dotted 
lines shift to the right. If at some redshift slabs of higher N HI had already become cold 
and unstable but the background flux intensity receives an extra input, the whole cooling 
process might stop. This can explain the lack of a large population of faint star forming 
disks of very low column density at z — 0. On the other hand, the observed deficit can also 
be due to a truncation of the distribution function of protodisks of low column densities. 

How does the theoretical scheme of Fig. 2 compare with observations? A precise answer 
to this question is difficult since the parameters N HI and Voo/Rq of protogalactic disks are 
unknown. However, as an illustration, we can use the results of de Blok et al. (1996) on 
Low Surface Brightness galaxies (LSBs) which are reminiscent of our models of LSPs. We 
have estimated Nhi and Voo/Rq for ten galaxies of their sample with slowly rising rotation 
curves, excluding LSBs with steeply declining gas column density profiles in the inner 
regions. The selected objects correspond to LSBs with 10 < -R25 < 30 arcsec and the results 
are shown by the filled triangles of Fig. 2. The data fall within the boundaries of the cold 
& unstable region and indicate that if LSPs were similar to LSBs they had the adequate 
parameters to follow the evolution described in this Letter. Within this region, systems of 
low column density and dynamical mass are those for which the ensuing starburst easily 
disrupts the gaseous central region and reduces considerably the gas mass content. Thus, 
further star formation is inhibited and the systems experience a fading away phase (CS). 
The identification of LSPs with the progenitors of LSBs is quite appealing since some of 
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the properties of the blue galaxies observed at moderate z can be explained if these are 
edge-on manifestations of LSBs (Dalcanton & Shectman 1996). Present-day LSBs might 
descend from those LSPs which, for environmental effects or some internal feature (higher 
dark matter or gas content), had a slower, and less violent star formation phase. 

We are grateful to E. E. Salpeter and to the referee for their comments and to R. 
Baglioni for his help with the figures. This work was supported in part by ASI-95-RS-120. 
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Table 1: List of reactions 

# Reaction 

1 H+ + e -> H + /«/ 

2 H + e -> H + hv 

3 H~ + H -> H 2 + e 

4 H- + H+ -> 2H 

5 H 2 + H+ -> + H 

6 H + hv -> H+ + e 

7 H- + /w/ -> H + e 

8 H 2 + /«/ -> + e 

9 H 2 + /«/ -> 2H 



Reference 

Seaton 1959 

De Jong 1972 

Bieniek 1980 

Peterson et al. 1971 

Hollenbach & McKee 1989 

Cen 1992 

Wishart 1979 

O'Neil & Reinhardt 1978 

see text 
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FIGURE CAPTIONS 

Fig. 1.— (a) Equilibrium gas temperatures as function of redshift for three cases: 
(1) N HI = 10 20 - 6 cm' 2 and A = B = 1, (2) N HI = 10 20 - 4 cm" 2 and A — B — 1, and 
(3) Nhi = 10 20 ' 4 cm" 2 and A = 10, B = 1 (continuos lines). The dashed lines show the 
time evolution of the temperature once the slab reaches the transition point marked with 
the asterisk. (6) Fractional abundances of H 2 after the transition point and prior to z = 
for the same cases as in (a). For all cases shown V^/Rq = 15 km s" 1 kpc" 1 . 

Fig. 2. — Regions in the N#/-l / 00 /i?o plane corresponding to the different thermal and 
dynamical regimes discussed in the text. The units of Voo/Rq are km s" 1 kpc" 1 , Nhi is in 
cm" 2 . The dotted lines connect points which become gravitationally unstable at the same 
redshift and are labelled accordingly. 
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